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The discovery of a low mass, pre-main-sequence stellar 
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ABSTRACT 

We report the serendipitous discovery of a population of low mass, pre-main sequence 
stars (PMS) in the direction of the Wolf-Rayet/O-star binary system j 2 Vel and 
the Vela OB2 association. We argue that j 2 Vel and the low mass stars are truly 
associated, are approximately coeval and that both are at distances between 360- 
490 pc, disagreeing at the 2a level with the recent Hipparcos parallax of 7 2 Vel, but 
consistent with older distance estimates. Our results clearly have implications for the 
physical parameters of the 7 2 Vel system, but also offer an exciting opportunity to 
investigate the influence of high mass stars on the mass function and circumstellar 
disc lifetimes of their lower mass PMS siblings. 

Key words: stars: pre-main-sequence - X-rays: stars - stars: Wolf-Rayet. 



1 INTRODUCTION 

As little as 15 years ago it was believed that high mass 
stars formed mainly in OB associations and that low mass 
stars formed mainly in T associations. In the last few years 
it has become increasingly recognized that the majority of 
low mass stars in the solar neighbourhood are in fact likely 
to have formed in OB associations. The mass function in 
OB associations may match popular log-normal parameter- 
isations of the general field population such as those pro- 
posed by Miller & Scalo (1979), with many low mass stars 
formed for every O/B star. The main reason for this shift in 
paradigm is the discovery of numerous low mass pre-main 
sequence (PMS) stars in young OB associations by virtue 
of their high levels of X-ray activity (see for example Wal- 
ter et al. 1994, 1999; Naylor & Fabian 1999, Preibisch & 
Zinnecker 1999 and references therein). Recently, some evi- 
dence has been found that these X-ray selected groups can 
be quite small, concentrated around just one or two high 
mass stars, which are themselves part of larger associations. 
Examples include /3 Cru (Feigelson & Lawson 1997), a Ori 
(Walter, Wolk & Sherry 1998) and r\ Cha (Mamajek, Lawson 
& Feigelson 1999). 

The investigation of these low mass PMS stars is of 
prime importance in establishing the form of the initial mass 
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function. Although OB associations are by definition un- 
bound, they are young enough that dynamical effects such 
as mass segregation and preferential evaporation will not 
have occurred; and they are usually free of the heavy, vari- 
able extinction that plagues observations of active star form- 
ing regions. If the majority of low mass stars are born in 
OB associations it is crucial to establish the influence that 
high mass neighbours have on the formation and evolution of 
their lower mass siblings. The winds and ionizing radiation 
of hot stars could influence the mass function and circum- 
stellar disc lifetimes of the lower mass stars, with implica- 
tions for angular momentum evolution and planet forma- 
tion. These ideas have gained currency with the discovery 
of evaporating discs around PMS stars in the Orion nebula 
(McCaughrean & O'Dell 1996) and theoretical studies show- 
ing that discs could be ionized and evaporated by the UV 
radiation fields of O stars (Johnstone, Hollenbach & Bally 
1998). In more extreme circumstances, nearby supernova ex- 
plosions have been invoked both as a means of terminating 
low mass star formation (Walter et al. 1994) or triggering it 
(Preibisch & Zinnecker 1999). 

In this paper we report the discovery, by X-ray selec- 
tion, of a low mass stellar population that seems likely to be 
associated with the nearest example of a Wolf-Rayet (WR) 
star, 7 2 Velorum (HD 68273, HIP 39953, WR11). Like about 
half of the ~ 200 galactic WR stars known, it is a binary 
system (WC8+08) with an orbital period of 78.5 days and 
a massive, interacting stellar wind. There is currently some 
controversy concerning the distance to 7 2 Vel, which impacts 
upon the deduced luminosities, masses and mass loss rates 
from the system. It is important to get these parameters 
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Figure 1. ROSAT PSPC X-ray image of the region around 7 2 
Vel. The greyscale is such that black represents > 1 photon per 
5x5 arcsec pixel and the brightest source in the centre of the image 
is 7 2 Vel. Another 108 significant X-ray sources have been found 
in this image. The solid outline shows the location of our optical 
CCD survey. 



right because, as the nearest WR, 7 Vel is an extreme test 
of stellar evolution models and calibrates the absolute mag- 
nitudes of WR stars. The Hipparcos parallax yields a dis- 
tance of 2581^ pc to 7 2 Vel (Schaerer, Schmutz & Grenon 
1997, van der Hucht et al. 1997), in marked contrast to pre- 
vious distance estimates which place it at 350-450 pc. The 
larger distance is in better agreement with the mean dis- 
tance to the Vela OB2 association (410±12pc), of which 7 2 
Vel is the most massive proper-motion member (de Zeeuw 
et al. 1999). 

The presence of a low mass stellar association in the 
extreme environment around 7 2 Velorum offers an excellent 
empirical test of the possible influence of winds and ionizing 
radiation on low mass stars and their discs. Additionally 
it gives us a chance to measure the distance to 7 2 Vel by 
matching PMS isochrones to the low mass stars. 



2 X-RAY OBSERVATIONS 

We suspected the presence of a low mass association around 
7 2 Vel from the large surrounding population of X-ray point 
sources seen in ROSAT images, taken for a programme to 
investigate its interacting stellar winds (see Willis, Schild 
& Stevens 1995 and Fig.l). The X-ray observations of the 
WR were retrieved from the R OS AT public archive and con- 
sisted of 10 Position Sensitive Proportional Counter (PSPC) 
datasets and 2 datasets taken with the High Resolution Im- 
ager. In this initial paper we discuss only the more sensitive 
PSPC results, which contain the vast majority of the X-ray 
point sources. 

We used the Starlink-distributed ASTERix data reduc- 
tion package in our X-ray analysis (Allan & Vallance 1995). 
The 10 PSPC datasets were sorted into 1° x 1° images, se- 



lecting pulse height channels 11 to 240 (approximately 0.1- 
2.4 keV photons) and excluding times with anomalously high 
background rates. The images were centered on 7 2 Vel (the 
brightest source in each dataset) and then summed. The re- 
sulting image has dimensions of 720 x 720 5 arcsecond pixels 
and an effective on-axis exposure time of 25.3 ks. We used 
the Point Source Searching (PSS - Allan 1992) algorithm 
to search for sources by Cash-statistic maximisation. By as- 
suming the background to be zero we obtained a preliminary 
list of 104 X-ray sources which were masked out of the im- 
age. The masked image was then patched and smoothed with 
a 75 arcsec FWHM gaussian to create a background map. 
We then executed the PSS algorithm again and found 109 
sources above a pseudo-gaussian significance level of 4.5a", 
which corresponds roughly to 1 spurious detection in the X- 
ray image. The positions of the X-ray sources were corrected 
for errors in the satellite aspect solution by comparing the 
optical and X-ray positions of 7 2 Vel and four other bright 
stars from the CDS SIMBAD database. We applied shifts 
of 4.5 arcsec in RA and 13.8 arcsec in Dec to the X-ray 
positions. 



3 OPTICAL PHOTOMETRY 

CCD photometry of the X-ray field was obtained on 8 Febru- 
ary 1999 with the 0.9-m telescope at the Cerro Tololo Inter- 
American Observatory. A Tek 2048x2048 CCD was used 
to give a 13.5x13.5 arcmin 2 field of view. Eight overlap- 
ping fields around 7 2 -Vel were surveyed in B VI with short 
(20,10,10s) and long (200,100,100s) exposures, together with 
five fields from Landolt (1992) to determine zeropoints, 
colour terms and extinction coefficients. Several standards 
with V — I c > 2.5 were observed. Figure 1 illustrates the 
location of the fields around 7 2 Vel. The external accuracy 
of our photometry was determined to be around 0.02 mag 
on the BVI C system. 

Photometry was performed for each of the eight fields 
using an optimal aperture algorithm (Naylor 1998; Totten 
et al. in preparation), and the results combined to give an 
optical catalogue of sources. Astrometry was performed by 
comparison with star positions in version A2.0 of the USNO 
catalogue (Monet 1998) and we estimate positional accura- 
cies of around 0.3 arcsec. The final catalogue contains 20617 
individual objects with their magnitudes, colours and posi- 
tions. Figure 2a shows the V versus V — I c colour- magnitude 
diagram (CMD). The location of a possible PMS low-mass 
star association is clearly visible above the bulk of the back- 
ground contamination. To this optical catalogue, we added 
the positions of bright (V < 11) stars found in this region 
from the SIMBAD database. 

To establish an appropriate cross-correlation radius to 
use between the X-ray and optical source lists we modelled 
the cumulative number of X-ray sources that were correlated 
with an optical source with V < 19 (see Jeffries, Thurston & 
Pye 1997 for details) . Assuming a uniform spread of optical 
sources, we determined that there were 77 correlations (from 
83 X-ray sources inside the CCD survey) within 10 arcsec of 
X-ray positions, that 66 of these would be true counterparts 
to X-ray sources, 11 would be spurious correlations and that 
the la X-ray error circle was 3.7 arcsecs. Figure 2b shows 75 
sources that have an optical counterpart within 10 arcsecs 
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Figure 2. (Top) CMD for the region around 7 2 Vel. The solid 
lines in this and the plot below are isochrones at 1, 5 and lOMyr 
calculated from the models of D'Antona & Mazzitelli (1997). The 
dashed box indicates the region from which PMS stars were se- 
lected for Fig. 3. (Bottom) The location of the X-ray source cor- 
relations in the CMD. Eleven spurious correlations are expected 
and are expected to lie in the clump at V ~ 19, V — J ~ 1. 
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Figure 3. (Top) The spatial distribution of PMS stars (from in- 
side the dashed box in Fig. 2) around 7 Vel. A clumping towards 
7 2 Vel (marked with a cross) is apparent. (Bottom) The relative 
spatial density of these PMS sources as a function of distance from 
7 2 Vel. This plot is normalized using a background population to 
account for non-uniform coverage in the CCD survey. 



(another two are bright stars without V — I colours). Clearly 
the X-ray emitting population coincides with the proposed 
PMS population in the CMD. Indeed, if we were to consider 
just a subset of the optical catalogue consisting of a broad 
strip containing all these PMS sources, we would only expect 
1 of these correlations to be spurious. 

We have calculated the X-ray properties of this pop- 
ulation and these along with the HRI observations will be 
reported in a subsequent publication. Briefly, the X-ray to 
bolometric flux ratio of these objects lies in the range 10~ 5 
to 10 -2 , broadly what we would expect from a population 
of young PMS stars. The cut-off in the PMS X-ray cor- 
relations at V ~ 18 is almost certainly due to the X-ray 
sensitivity. To be detected in X-rays, fainter objects would 



have to have higher than feasible X-ray to bolometric flux 
ratios. However, it is clear that the PMS sequence we have 
found extends down to the limits of our optical survey at 
V ~ 20.5. 



4 DISCUSSION 

The appearance of Fig. 2 should leave the reader in no doubt 
that we have found a young and exceptionally rich popula- 
tion of low mass active stars in the direction of 7 2 Vel. The 
central question to be answered is whether these sources 
are physically close to 7 2 Vel and/or whether they are back- 
ground members of the Vela OB2 association. We can tackle 
this problem in a number of ways. 
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4.1 Are the PMS stars and j 2 Vel aligned? 

Figure 3a shows the spatial distribution of PMS stars se- 
lected from the CMD in a strip enclosing the bulk of the 
X-ray sources (marked with a dashed box in Fig.2). This 
box was chosen to avoid background contamination. We de- 
termined the radial distribution of these stars, centred on 7 2 
Vel. The distribution is normalized using the radial distribu- 
tion of background stars with a similar V magnitude range, 
but a colour range of 0.8< V — I c <1.7, under the assump- 
tion that the background stars are uniformly distributed. 
The resulting radial distribution is shown in Fig. 3b, which 
exhibits a small but significant (at the 3<r level) peak within 
5 arcmin of the centre of the CCD survey and -y 2 Vel. The 
point closest to 7 2 Vel is missing because the 30 arcsec re- 
gion immediately surrounding 7 2 Vel is swamped by its light 
and no accurate photometry was obtained there. In a sim- 
ilar fashion we can show that the X-ray sources are also 
marginally concentrated toward the centre of the field even 
after correction for the PSPC vignetting function. 

There is thus some evidence that the PMS stars and 7 2 
Vel are spatially correlated, although this does not rule out 
a chance alignment of 7 Vel with a background cluster of 
low mass stars in Vela OB2. In particular, it is possible that 
any concentration we see could be associated with 7 1 Vel, a 
B2III, single lined spectroscopic binary which is a common 
proper motion companion to 7 2 Vel. 7 1 Vel is also a likely 
member of the Vela OB2 association, has a spectroscopic 
parallax of ~ 450 pc and is located only 41 arcsecs at a PA 
of 220° from 7 2 Vel. 

4.2 Are the PMS stars and 7 2 Vel at the same 
distance? 

The isochrones plotted in Fig.2 come from the D'Antona & 
Mazzitelli (1997) low mass evolutionary models. We have 
converted from bolometric luminosity and effective temper- 
ature to magnitude and colour using empirical bolometric 
corrections as a function of colour and a colour-effective 
temperature relationship derived by forcing the well stud- 
ied low mass stars in the Pleiades cluster to fit a 125 Myr 
isochrone (see Jeffries & Tolley 1998). The assumed distance 
and reddening are those appropriate for the Vela OB2 asso- 
ciation of 410 pc and E(V — I) ~ 0.06. It is clear that if the 
PMS stars are at the mean distance of the Vela OB2 asso- 
ciation they appear to be about 4 ± 2 Myr old (taking into 
account the likely presence of unresolved binary systems). 
If however, the PMS association were at the Hipparcos dis- 
tance to 7 2 Vel, the isochrones would be shifted upwards by 
1 magnitude. In that case we would deduce that the PMS 
population had an age > 10 Myr. The age of 7 2 Vel, based 
upon the mass of the O star (~ 3OM0) at the Hipparcos dis- 
tance, is less than 5 Myr (Schaerer et al. 1997; de Marco & 
Schmutz 1999). Thus if the Hipparcos distance is adopted, 
the PMS stars and 7 2 Vel cannot be at the same distance 
and coeval. However, if 7 2 Vel were at the mean distance 
of the Vela OB2 association, it would be more massive (see 
below), slightly younger and could easily be coeval with the 
PMS population. An age range of 2-6 Myr would be com- 
patible with distances between 490 and 360 pc. 



4.3 How far away is 7 2 Vel? 

The evidence that 7 2 Vel is as close as 258 pc from the 
Hipparcos parallax should be treated with some caution, 
de Zeeuw et al. (1999) examine the high mass membership 
of the Vela OB2 association on the basis of both proper 
motions and parallaxes. The association is well defined by 
proper motions and has an angular radius of about 6 de- 
grees. The mean parallax from 93 members corresponds to 
410 ± 12 pc. The parallax dispersion can be quite well mod- 
elled as a gaussian with a a — 0.68 mas, such that ~ 95% of 
the members appear to lie between 260 pc and 900 pc. The 
dispersion is reasonably consistent with the errors on the 
individual points and not inconsistent with the idea that all 
the stars are at nearly the same distance. Indeed, if we were 
to assume that the front to back size of Vela OB2 were sim- 
ilar to its diameter on the sky, then all the stars should be 
contained within ±40 pc. We would therefore interpret the 
parallax to 7 2 Vel simply as a ~ 2a deviation and that its 
distance was 410 ± 40 pc. 

The companionship of 7 1 Vel is also evidence for a dis- 
tance closer to the mean Vela OB2 distance. Although these 
objects have a common proper motion, they were too close 
together on the sky for Hipparcos to obtain independent 
parallaxes. The distance to 7 1 Vel from its spectral type 
and photometry is almost certainly 400-500 pc (e.g. Abt et 
al. 1976; Hernandez & Sahade 1980). The idea that 7 2 Vel 
is a foreground object randomly placed within 1 arcminutc 
of another bright member of the Vela OB2 association has 
a probability of only ~ 10 -3 . There is also reasonable evi- 
dence that the radial velocities of the two systems are sim- 
ilar. Hernandez & Sahade (1980) quote 9.7 ± l.Okms -1 for 

7 1 Vel and Niemela & Sahade (1980) give 12 ± 1 kms" 1 for 

7 2 Vel, although it is likely that the accuracy of the latter 
result is exaggerated (Schmutz et al. 1997). 

4.4 Can the PMS stars be at a range of distances? 

Preibisch & Zinnecker (1999) observed a very similar PMS 
CMD in the Upper Sco OB association, with a similar ver- 
tical scatter of about ±0.6 mag about the isochrones. They 
showed that if one takes into account unresolved binaries, 
photometric errors and allowed a ~ 10% range in distance, 
that the scatter around the PMS isochrones was exactly as 
expected for a coeval population at ~ 5 Myr. We therefore 
similarly conclude that the CMD in Fig. 2 shows no evidence 
for a large spread in either the age or distance of the PMS 
population we have found. In particular, we can rule out 
any distance modulus spread in a coeval population that is 
larger than a few tenths of a magnitude, or any age spread 
in a co-spatial population of more than a Myr or so (for a 
mean age of ~ 4 Myr) . Thus unless there is a conspiracy to 
place older stars closer to us, the PMS association seems 
likely to have a relatively narrow spread around an age and 
distance that are incompatible with the deduced age for 7 2 
Vel and its Hipparcos distance. 

4.5 Is a larger distance to 7 2 Vel consistent with 
its physical properties? 

Our findings challenge the conclusions of several recent pa- 
pers which use the Hipparcos parallax of 7 2 Vel and its er- 
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ror to derive: the absolute magnitude of the system and its 
components; system masses from the interferometric binary 
separation of Hanbury- Brown (1970); the O star luminosity, 
mass and age from stellar evolution models and hence the 
orbital inclination and further mass estimates from radial 
velocity curves (see van der Hucht et al. 1997; Schaerer et 
al. 1997, Schmutz et al. 1997, de Marco & Schmutz 1999). 
A distance as large as 410 pc for j 2 Vel significantly changes 
the system parameters deduced in these papers. The system 
luminosity increases by a factor 2.5. The effective tempera- 
ture and luminosity deduced for the O star would then give 
it a mass > 40 Mq and an age < 3Myr, compared with the 
values of 30 Mq and 3.6 Myr quoted by de Marco & Schmutz 
(1999). At this larger distance, an age of 2-3 Myr could be 
compatible with the low mass PMS stars we have found. 
The absolute magnitude of the O star would decrease to 
—6.0 ± 0.3 (van der Hucht et al. 1997), which argues for a 
supergiant rather than a giant classification. Van der Hucht 
et al. comment that this would be in better agreement with 
published spectra of Conti & Smith (1972) and Niemela & 
Sahade (1980). As the mass ratio from the radial velocity 
curves is fixed, the WR mass increases by a similar frac- 
tion. The total system mass, based on a binary separation 
of 4.3 ± 0.5 mas, increases from 30 ± 10 M to 120 ± 40 M 
(now comfortably exceeding the minimum mass from radial 
velocity curves - Niemela & Sahade 1980) and the binary 
inclination is reduced to around 50° to explain the radial 
velocity curves. 



5 CONCLUSIONS 

From our discussion there seem to be two possible scenar- 
ios. (1) That the PMS stars are approximately at the same 
distance and age as j 2 Vel, and that this distance places 7 2 
Vel within the Vela OB2 association at 360-490 pc. (2) That 
the PMS stars are part of the Vela OB2 association, possibly 
surrounding 7 1 Vel, but that 7 2 Vel is an isolated foreground 
object with no surrounding low mass stars at a similar age. 
We believe that (1) is far more plausible than (2) because 
of the dispersion in the Vela OB2 Hipparcos parallaxes and 
the likely association of 7 1 and 7 2 Vel. Recently, the idea 
that 7 2 Vel could form in isolation without accompanying 
low mass stars has also been challenged by the near IR de- 
tection of a K-type PMS companion only 4.7 arcsec distant 
(Tokovinin et al. 1999). 

If the low mass PMS stars we have found are truly in 
the vicinity of 7 2 Vel, they represent an exciting opportu- 
nity to explore the influence of adjacent high mass loss stars 
and ionizing UV radiation fields on the mass function and 
circumstellar disc lifetimes of low mass stars. It will be in- 
teresting to compare the frequencies of T-Tauri discs around 
these stars with the frequencies found in T associations and 
OB associations with similar ages. The PMS stars in Fig. 2 
have masses, found from the D'Antona & Mazzitelli (1997) 
models, down to (an age dependent) mass of about 0.15 Mq. 
The mass function will be addressed when we have a better 
census of the association membership. 
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